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Abstract 

We present next-to-leading order predictions for double transverse-spin asymmetries 
in Drell-Yan dilepton production initiated by proton-antiproton scattering. The 
kinematic region of the proposed PAX experiment at GSI: 30 < s < 200 GeV 2 
and 2 < M < 7 GeV is examined. The Drell-Yan asymmetries turn out to be 
large, in the range 20-40%. Measuring these asymmetries would provide the cleanest 
determination of the quark transversity distributions. 
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1. The experiments with antiproton beams planned for the next decade in 
the High-Energy Storage Ring at GSI will provide a variety of perturbative 
and non-perturbative tests of QCD [1], In particular, the possible availability 
of transversely polarised, antiprotons opens the way to direct investigation 
of transversity, which is currently one of the main goals of high-energy spin 
physics [2], The quark transversity (i.e. transverse polarisation) distributions 
Aqrq were first introduced and studied in the context of transversely polarised 
Drell Yan (DY) production [3]; this is indeed the cleanest process probing 
these quantities. In fact, whereas in semi-inclusive deep-inelastic scattering 
transversity couples to another unknown quantity, the Collins fragmentation 
function [4], rendering the extraction of A Y q a not straightforward task, the 
DY double-spin asymmetry 

dy = dff rT - ckx n A to_ 

TT ~ dull - dun aunp ( j 

only contains combinations of transversity distributions. At leading order, for 
instance, for the process p^p 1 ' —> £ + £~X one has 

DY _ Eg e*[A T q(xi, M 2 ) A T q{x 2 , M 2 ) + A T q(xi, M 2 ) A T q{x 2 , M 2 )] 

Eg e\[q{x u M 2 ) q(x 2 , M 2 ) + q(x 1 , M 2 ) q(x 2 , M 2 )} 

(2) 

where M is the invariant mass of the lepton pair, q(x, M 2 ) is the unpolarised 
distribution function, and ott is the spin asymmetry of the QED elementary 
process qq —> t + l~. In the dilepton centre-of-mass frame, integrating over the 
production angle 6 , one has 


a T r (</?) = \ cos 2(p, (3) 

where ip is the angle between the dilepton direction and the plane defined by 
the collision and polarisation axes. 

Measurement of ptpt DY is planned at RHIC [5]. It turns out, however, that 
Att(pp) is rather small at such energies [6,7,8], no more than a few percent 
(similar values are found for double transverse-spin asymmetries in prompt- 
photon production [9] and single-inclusive hadron production [10]). The rea¬ 
son is twofold: 1) (pp) depends on antiquark transversity distributions, 
which are most likely to be smaller than valence transversity distributions; 2) 
RHIC kinematics (i /s = 200 GeV, M < 10 GeV and X\X 2 = M 2 / s < 3 x 10 -3 ) 
probes the low-rc region, where QCD evolution suppresses A Y q(x, M 2 ) as com¬ 
pared to the unpolarised distribution q(x, M' 2 ) [11,12], The problem may be 
circumvented by studying transversely polarised proton-antiproton DY pro¬ 
duction at more moderate energies. In this case a much larger asymmetry 
is expected [6,13,14] since (pp) is dominated by valence distributions at 
medium x. The PAX collaboration has proposed the study of jdp Drell-Yan 
production in the High-Energy Storage Ring (HESR) at GSI, in the kinematic 
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region 30 GeV 2 < s < 200 GeV 2 , 2 GeV < M < 10 GeV and xix 2 > 0.1 [15]. 
An antiproton polariser for the PAX experiment is currently under study [16]: 
the aim is to achieve a polarisation of 30-40%, which would render the mea¬ 
surement of Ayr (pp) very promising. 

Leading-order predictions for the pp asymmetry at moderate s were presented 
in [13]. It was also suggested there to access transversity in the J/tp resonance 
production region, where the production rate is much higher. The purpose of 
this paper is to extend the calculations of [13] to next-to-leading order (NLO) 
in QCD. 1 This is a necessary check of the previous conclusions, given the 
moderate values of s in which we are interested. We shall see that the NLO 
corrections are actually rather small and double transverse-spin asymmetries 
are confirmed to be of order 20-40%. 


2. The kinematic variables describing the Drell-Yan process are (1 and 2 
denote the colliding hadrons): 


fi = VTe 2/ , 6 = y = ^ In ^, (4) 

with r = M 2 / s. We denote by X\ and x 2 the longitudinal momentum fractions 
of the incident partons. At leading order, 0 and £ 2 coincide with X\ and x 2 , 
respectively. The QCD factorisation formula for the transversely polarised 
cross-section for the proton-antiproton Drell-Yan process is 


(IAtct 

d M d y dp 



dx 1 dx 2 [A T q(x 1 ,p 2 )A T q(x 2 ,p 2 ) 


+ Axq(xi, p 2 )A T q(x 2 , p 2 ) 


dA r d 


d M dy dp ’ 


( 5 ) 


where // is the factorisation scale and we take the quark (antiquark) distri¬ 
butions of the antiproton equal to the antiquark (quark) distributions of the 
proton. Note that, since gluons cannot be transversely polarised (there is no 
such thing as a gluon transversity distribution for a spin one-half object like the 
proton), only quark-antiquark annihilation subprocesses (with their radiative 
corrections) contribute to dA^cr. In Eq. (5) we use the fact that antiquark dis¬ 
tributions in antiprotons equal quark distributions in protons, and viceversa. 
At NLO, i.e. at order a s , the hard-scattering cross-section dAxA^, taking 


1 The results presented here were communicated at the QCD-PAC meeting at GSI 
(March 2005) and reported by one of us (M.G.) at the Int. Workshop “Transversity 
2005” (Como, September 2005) [17]. 
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the diagrams of Fig. 1 into account, is given by [7] 


dATct (1 ^ MS 2a 2 a s (/i 2 ) 4t(xiX2 + r) 

dM dy d<^ 9 sM F 


x S <S(a?i - ^i)A(a: 2 - 6) 

1 


2tt XiX 2 (Xi +6) (^2 + 6) 

2(l-6)(l-6) , vr : 


cos(2<^) 


+ <5(^1 — ■Ci) 


. In 

4 r 

2x 2 (1-6) 


+ T “ 2 
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ln(x 2 - 6)' 
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(® 2 - 6 )+ t(x 2 + ^ 2 ) V 
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+ 


1 1 ^ 

t i b-T ln — 

x 2 - K2 J + x 2 — £ 2 x 2 
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(®2 - 6)- 


+ 


( 6 ) 


where we have taken the factorisation scale // equal to the renormalisation 
scale. In our calculations we set /i = M. 



Fig. 1. Elementary processes contributing to the transverse Drell-Yan cross-section 
at NLO: (a, b) virtual-gluon corrections and (c) real-gluon emission. 

The unpolarised Drell-Yan differential cross-section can be found, for instance, 
in [18]; besides the diagrams of Fig. 1, it also includes the contribution of 
quark-gluon scattering processes. 


3. To compute the Drell-Yan asymmetries we need an assumption for the 
transversity distributions, which as yet are completely unknown. We might 
suppose, for instance, that transversity equals hclicity at some low scale, 
as suggested by confinement models [11] (this is exactly true in the non- 
relativistic limit). Thus, one possibility is 

A T q(x,nl) = Aq(x,fi 2 0 ) , ( 7 ) 

where typically /r 0 < 1 GeV. Another possible assumption for A xq is the 
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saturation of Soffer’s inequality [19], namely 


A T q(x,nl) 


\ q{x,nl) + kq{x,nl) , 


which represents an upper bound on the transversity distributions. 


( 8 ) 


Since Eqs. (7) and (8) make sense only at very low scales, in practical calcula¬ 
tions one has to resort to radiatively generated helic-ity and number densities, 
such as those provided by the GRV fits [20]. The GRV starting scale is indeed 
(at NLO) /i ( j = 0.40 GeV 2 . We should however bear in mind that in the GRV 
parametrisation there is a sizeable gluon contribution to the nucleon’s helic-ity 
already at the input scale (A g is of order 0.5). On the other hand, as already 
mentioned, gluons do not contribute to the nucleon’s transversity. Thus, use 
of Eq. (7) with the GRV parametrisation may lead to substantially underes¬ 
timating the quark transversity distributions and hence is a sort of “minimal 
bound” for transversity. Incidentally, the experimental verification or other¬ 
wise of the theoretical predictions of A TT based on the low-scale constraints 
(7, 8) would represent an indirect test of the “valence glue” hypothesis behind 
the GRV fits. Note too that, although the assumption (7) may, in principle, 
violate the Soffer inequality, we have explicitly checked that this is not the 
case with all the distributions we use. 


After setting the initial condition (7) or (8), all distributions are evolved at 
NLO according to the appropriate DGLAP equations (for transversity, see [21]; 
the numerical codes we use to solve the DGLAP equations are described 
in [22]). The u sector of transversity is displayed in Fig. 2 for the minimal 
bound (7) and for the Soffer bound (8). 

The transverse Drell-Yan asymmetry integrated over M between 

2 GeV and 3 GeV (i.e. below the J/ip resonance region), for various values 
of s is shown in Fig. 3. As can be seen, the asymmetry is of order of 30% 
for s = 30 GeV 2 (fixed-target option) and decreases by a factor two for a 
centre-of-mass energy typical of the collider mode (s = 200 GeV 2 ). The cor¬ 
responding asymmetry obtained by saturating the Soffer bound, that is by 
using Eq. (8) for the input distributions, is displayed in Fig. 4. As expected, 
it is systematically larger, rising to over 50% for fixed-target kinematics. 

Above the J/ip peak A ^ /(Itt appears as shown in Fig. 5, where we present 
the results obtained with the minimal bound (7). Comparing Figs. 3 and 5, we 
see that the asymmetry increases at larger M (recall though that the cross- 
section falls rapidly with growing M ). 

The importance of NLO QCD corrections may be appreciated from Fig. 6, 
where one sees that the NLO effects hardly modify the asymmetry since the K 
factors of the transversely polarised and unpolarised cross-sections are similar 
to each other and therefore cancel out in the ratio. As for the dependence on 
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Fig. 2. The u and u transversity distributions, as obtained from the GRV parametri- 
sation and Eq. (7), top panel, or Eq. (8), bottom: xA tu at /a 2 = /j,q = 0.40GeV 2 
(dashed curve) and /a 2 = 9GeV 2 (solid curve); xA t'U at n 2 = Hq = 0.40 GeV 2 
(dotted curve) and fi 2 = 9GeV 2 (dot-dashed curve). Note that the u transversity 
distributions have been multiplied by a factor of 10. 

the factorisation scale /j (we recall that the results presented in all figures are 
obtained setting // = M), we have repeated the calculations with two other 
choices (/i = 2 M and /i = M/2 ) and found no sensible differences. 

A caveat is in order at this point. The GSI kinematics is dominated by the do¬ 
main of large r and large z = r/x 1 X 2 , where real-gluon emission is suppressed 
and where there are powers of large logarithms of the form ln(l — z), which 
need to be resummed to all orders in a s [23]. It turns out that the effects 
of threshold resummation on the asymmetry A^ in the regime we are con¬ 
sidering, although not irrelevant, are rather small (about 10%) if somewhat 
dependent on the infrared cutoff for soft-gluon emission. 
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Fig. 3. The NLO double transverse-spin asymmetry Att(d)/(Itt, integrated be¬ 
tween M = 2 GeV and M = 3 GeV, for various values of s; the minimal bound (7) 
is used for the input distributions. 



y 

Fig. 4. As Fig. 3, but with input distributions corresponding to the Soffer bound (8). 

The feasibility of the A TT measurement at GSI has been thoroughly investi¬ 
gated by the PAX Collaboration (see App. F of [15]). In collider mode, with 
a luminosity of 5 • 10 30 cm _2 s _1 , a proton polarisation of 80%, an antipro¬ 
ton polarisation of 30% and considering dimuon invariant masses down to 
M — 2 GeV, after one year’s data taking one expects a few hundred events 
per day and a statistical accuracy on A TT of 10-20%. 

4. Before concluding, we briefly comment on the possibility of accessing 
transversity via .//'</; production in pp scattering. It is known that the dilep¬ 
ton production rate around M — 3 GeV, i.e. at the J/p> peak, is two orders of 
magnitude higher than in the region M ~ 4 GeV. Thus, with a luminosity of 
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Fig. 5. The NLO double transverse-spin asymmetry Att(d)/o,tt, integrated be¬ 
tween M = 4GeV and M = 7 GeV, for various values of s; the minimal bound (7) 
is used for the input distributions. 
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Fig. 6. NLO vs. LO double transverse-spin asymmetry Att{u)/&tt at M = 4 GeV 
for s = 45 GeV 2 and s = 200 GeV 2 ; the minimal bound (7) is used for the input 
distributions. 

5 • 10 30 cm _2 s _1 , one expects a number of pp —> J/tf} —> £~£ + events of order 
10 5 per year at GSI collider energies. This renders the measurement of Att 
in the J/ip- resonance region extremely advantageous from a statistical point 
of view. 

As explained in [13], if J/'ijj formation is dominated by the qq annihilation 
channel, at leading order the double transverse-spin asymmetry at the J/ijj 
peak has the same structure as the asymmetry for Drcll-Yan continuum pro¬ 
duction, since the J/ip is a vector particle and the qq J coupling has the 
same hclicity structure as the qq'y* coupling. The CERN SPS data [24] show 

















































































































































































































































that the pp cross-section for J/i/j production at s = 80 GeV 2 is about ten 
times larger than the corresponding pp cross-section, which is a strong indi¬ 
cation that the (/(/-fusion mechanism is indeed dominant. Therefore, at the s 
values of interest here (s < 200 GeV 2 ) dilepton production in the J/pj reso¬ 
nance region can be described in a manner analogous to Drcll-Yan continuum 
production, with the elementary subprocess qq —> 7 * —> £~£ + replaced by 
qq —> J/if) —> £~£ + [25]. Using this model, which successfully accounts for the 
SPS J/ip production data at moderate values of s, it was found in [13] that 
the transverse asymmetry at the J/^ peak is of the order of 25-30%. 


At next-to-leading order, due to QCD radiative corrections, one cannot use 
a point-like qq J/ip coupling, and therefore it is not possible to extend in a 
straightforward way the model used to evaluate Att at leading order. Were 
NLO effects not dominant, as is the case for continuum production, one could 
still expect the J/ip asymmetry to be quite sizeable, but this is no more than 
an educated guess. What we wish to emphasise, however, is the importance 
of experimentally investigating the J/ip double transverse-spin asymmetry, 
which can shed light both on the transversity content of the nucleon and on 
the mechanism of J/ip formation (since gluon-initiated hard processes do not 
contribute to the transversely polarised scattering, the study of Att in the 
J/ip resonance region may give information on the relative weight of gluon 
and quark-antiquark subprocesses in J/ip production). 


5. In conclusion, experiments with polarised antiprotons at GSI will repre¬ 
sent a unique opportunity to investigate the transverse polarisation structure 
of hadrons. The present paper, which confirms the results of [13], shows that 
the double transverse-spin asymmetries are large enough to be experimentally 
measured and therefore represent the most promising observables to directly 
access the quark transversity distributions. 
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